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PHASE-DIVERSITY WAVEFRONT SENSOR 



The present invention relates to a novel wavefront sensor having a novel aberration 
means which can be constructed as a diffractive optical element (DOE) and which is 
5 suitable for use in adaptive optics. 

Adaptive optics (AO) is a technique in which real-time modification of the properties 
of an optical system based on measured optical output is used to indicate system 
optical performance. AO systems are widely used for compensation of turbulence in 
10 imaging through turbulent media in astronomical 1 and military 11 applications, for 
control of laser-beam delivery 111 , for control of laser resonator properties lv , for 
ophthalmic applications v and in other areas. In addition, aberration means can be used 
for measuring the distorting effect of a mirrored surface on an incident wavefront 
when it is reflected from the mirrored surface. 

15 

AO systems typically comprise 3 components: 

• a wavefront modulator (WFM) which alters the optical properties of the system in 
response to a command signal; 

• a wavefront sensor (WFS) which monitors the difference in state between the 
20 desired optical performance and the current optical performance of the system; 

and 

• a control loop which drives the WFM in response to the output from the (WFS). 



For control of an adaptive optical system it is not necessary to reconstruct the input 
25 wavefront. A sufficient condition for satisfactory operation of an adaptive optical 
system is the ability to drive a wavefront modulator using a null sensor. In these 
circumstances, a control signal derived from a wavefront sensor system is used to 
indicate the size and, preferably, the location and the direction of the wavefront error. 
Thus, if the wavefront modulator is providing foil correction of the input wavefront 
30 error, the control signal will be zero and the wavefront modulator will not be driven 
from its present position. The simplest example of such systems are based on multi- 
dither techniques applied, for example vl , to the correction of thermal lensing in laser 
resonators. 
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The phase-diversity™ algorithm can also be used to reconstruct the wavefront phase 
by using data corresponding to images of the input wavefront intensity on two planes 
normal to the direction of propagation and located at different positions along the axis 

5 of propagation. The approach taken is generally close to the c two-defocus' method 
used in microscopy™ and operates with two close-to-focus images* Although the 
data planes are generally described as symmetrically placed about the image plane, 
they can equally well be symmetrically placed about the system input pupil*, in which 
case the phase diversity algorithm becomes essentially the same as the wavefront 

10 curvature algorithm 5 ". 

For measurements on planes symmetrically spaced either side of either the image or 
the pupil planes, the intensity on the two measurement planes will be identical, and 
the difference between the images will be zero, if and only if, the wavefront in the 
15 entrance pupil plane is a plane wavefront. 

The measurement of the difference between the intensity on the two data planes (the 
phase-diverse data), thus satisfies the requirements for a null sensor because it 
provides a zero result when a plane wavefront occurs at the (input) pupil plane . If 
20 the input wavefront is distorted, the propagation between the measurement planes 
results in excess convergence (concave wavefront) or excess divergence (convex 
wavefront) and the resulting intensity difference between the measurement planes is 
indicative of the location, magnitude and direction of the wavefront curvature. 

25 Figure la shows a concave wavefront 3 propagating from left to right. For the 
sampled rays 3 shown, there is a local increase of intensity on one plane 5 and a 
decrease on the other plane 9. A convex wavefront 7 leads to a decrease in the local 
intensity. In each case, the intensities recorded in the +1 and -1 diffraction orders 
represent the intensities on the ends of the cylinder. 

30 

Figure lb shows that ambiguity can arise if the curvature of the wavefront 1 1 is so 
severe that a focus point occurs within the volume sampled by the measurements. 
The upper part of each figure leads to the same intensity increase. 
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Data on two image planes can be recorded using various approaches, including 
physical displacement of the image plane x , use of a vibrating spherically-distorted 
mirror, beam splitters and folded optical paths, or by the use of off-axis Fresnel 
lenses ,x . Each of these approaches allows the recovery of phase information from two 
5 data sets recorded under different focus conditions. As such, all of these approaches 
are related to the two-defocus methods applied in electron microscopy in the 1970s. 

Alternatively, a matched filter approach xu , can be adopted by using a diffractive 
optical element to provide a pair of spots whose axial intensity (passing through a 
10 pinhole) defines the amplitude of each aberration mode in a distorted wavefront. 

For reconstruction of the wavefront phase the inverse problem may be solved 
iteratively xm,xlv or presented in terms of the differential Intensity Transport Equation 
(ITE) and solved through the use of Green's functions™ with the phase-diverse data 
15 providing an estimate of the axial derivative of the intensity. Assumptions imposed 
through use of the ITE involve uniformity of the input intensity in the entrance pupil, 
continuity in the wavefront phase and continuity of the first derivative of the 
wavefront phase™. 

20 At the expense of assumptions imposed on the uniformity of the input illumination, 
the work on phase-diversity has progressed a long way from the early electron 
microscope applications of the two-defocus technique and ITE-based approaches have 
demonstrated real-time data reduction with high (sub-nanometre) accuracy™ 1 . 
However, in all cases the two data sets are recorded under conditions where the 

25 wavefront is subject to a known defocus aberration between the two measurements. 

The restrictions on the uniformity and continuity of the input wavefront mean that 
current wavefront sensors cannot be used with scintillated, discontinuous, or multiply- 
connected wavefronts. There are many areas in which a generalised method which can 
30 cope with scintillated, discontinuous and disconnected wavefront segments would be 
a great advantage. One example is in polishing machines, where it would be useful to 
measure in situ the surface profile of the object being polished. This would inevitably 
involve laser illumination of rough surfaces. It would also be useful to be able to 
image silicon circuitry, which by definition, is highly discontinuous. In addition. 
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modern large-aperture astronomical telescopes are likely to be constructed using 
segmented, thus multiply-connected, mirror segments. 

It is an object of the present invention to provide a wavefront sensor which uses the 
5 null sensor concept that will produce a null output if the input wavefront is a plane 
wave, but which will return an error signal when the input is distorted. A plane wave 
is defined as having constant phase across its wavefront. Without loss of generality 
that phase can be taken to be zero, so that its Fourier transform of the input 
wavefunction is Hermitian. 

10 

In accordance with a first aspect of the present invention there is provided a 
measuring apparatus for determining data relating to the shape of an input radiation 
wavefront, the wavefront shape being describable at a pre-determined location in an 
optical system: 

15 the apparatus comprising aberration means, the shape of which is defined by a filter 
function; 

detection means having a radiation sensitive surface capable of detecting the intensity 
of incident radiation on the surface, the detection means being coupled to an output 
device that provides a measure of the intensity of the incident radiation; 
20 wherein the aberration means acts on any input wavefront shape to produce first and 
second output radiation signals that in combination provide data from the output 
device on the extent to which the wavefront shape is non-planar. 

Preferably, the aberration means creates at least two filter functions, these filter 
25 function being complex conjugate pairs. 

Preferably, the filter functions are complex valued and have non-mixed symmetry. 
This feature allows the invention to be used for the analysis of scintillated 
discontinuous or disconnected wavefront segments. 

30 

The filter function cannot be described by a pure defocus or an aberration means 
consisting of a diffractive optical element with essentially quadratic curvature in the 
grating. Defocus diffractive optical elements cannot be used in the analysis of 
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scintillated, discontinuous or disconnected wavefront segments without generating 
boundary errors related to the wave shape. 

Preferably, the measuring apparatus further comprising a wavefront modulator. 

5 

Preferably, the wavefront modulator is configured to transform a wavefront being 
describable by means of a complex function into a wavefront being describable by a 
real function. 

10 Preferably, the output device is provided with calculation means for calculating the 
difference between the first and second radiation signals. 

Preferably, the radiation sensitive surface of the detection means is provided with 
elements that allow the measurement of radiation intensity at different points across 
15 the surface of the detection means. 

Preferably, the output device is provided with calculation means for calculating the 
difference between the first and second radiation signals at different points across the 
surface of the detection means. 

20 

Preferably, the wavefront modulator is coupled to the output device such that the 
wavefront modulator is distorted to provide a correction to a non-planar input 
radiation wavefront. 

25 Preferably, the aberration function is a weighted sum of Zernike polynomials arranged 
to equalise the signal generated from each mode of deformation in the input wavefront 
according to the expected statistical distribution of such modes in the input wavefront. 

Preferably, the aberration means is arranged such that the complex conjugate filter 
30 functions of the aberration means are associated with diffraction orders of the same 
order but having different sign. 



Preferably, the first and second output radiation are produced simultaneously. 
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Optionally, the first and second output radiation are produced sequentially. 

Preferably, the aberration means is a diffractive optical element. 

5 Optionally, the aberration means can be a variable-shape optical mirror 

Optionally, the aberration means is a variable refractive index device. 

Optionally, the variable refractive index device is a liquid crystal phase modulator 
10 used sequentially to provide complex conjugate aberrations or to provide a 
programmable diffractive optical element. 

Optionally, the aberration means is a deformed reflective surface where the 
illumination of that surface from each side produces the complex conjugate aberration 
15 functions. 

In accordance with a second aspect of the present invention there is provided a 
method for determining data relating to the shape of an input radiation wavefront, the 
wavefront shape being describable at a pre-determined location in an optical system: 
20 the method comprising the steps of 

transmitting said input radiation wavefront through an aberration means, the shape of 

which is defined by a filter function 

detecting the intensity of incident radiation on a surface, 

sending the detected intensity to an output device that provides a measure of the 
25 intensity of the incident radiation; 

wherein the aberration means acts on any input wavefront shape to produce first and 
second output radiation signals that in combination provide data from the output 
device on the extent to which the wavefront shape is non-planar. 

30 Preferably, the iterative reconstruction is used to reconstruct the wavefront shape. 



Preferably, the filter function is complex valued and has non-mixed symmetry. 
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In accordance with a third aspect of the present invention there is provided a 
diffractive optical element comprising aberration means, the shape of which is defined 
by a filter function that is complex valued and which has non-mixed symmetry. 

5 Preferably, the aberration means can be a variable-shape optical mirror 

Preferably, the aberration means is a variable refractive index device. 

More preferably, the variable refractive index device is a liquid crystal phase 
10 modulator used sequentially to provide complex conjugate aberrations. 

Preferably, the aberration means is a deformed reflective surface where the 
illumination of that surface from each side produces the complex conjugate aberration 
functions. 

15 

In accordance with a fourth aspect of the present invention there is provided a method 
for creating a diffractive optical element comprising the steps of defining aberration 
means the shape of which is defined by a filter function and which has non-mixed 
symmetry, applying said filter function to an optical element to create a diffractive 
20 optical element. 

Sufficient conditions for the null sensor are that it should give null output for input 
wavefunctions whose Fourier transforms have Hermitian symmetry, and that it give 
non-zero output for non-plane input wavefronts. In addition, the filter function should 
25 be complex. Hermitian symmetry exists where a function reflected in the origin is the 
complex conjugate of itself. 

Conditions where the real or imaginary part of the filter function is zero should be 
avoided to prevent a situation arising where the sensor will give null output for any 
30 input. In addition, mixed symmetries of the real and imaginary parts of the filter 

function should be avoided, since such filter functions provide an output signal even 
for a plane wave input 
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For unambiguous sensor output the real and imaginary parts of the filter function 
should have either both odd or both even symmetry. From these conditions it is 
possible to construct a null sensor without having placed any limiting assumptions on 
the input wavefront. This method will therefore be suitable for scintillated and 
5 discontinuous input wavefronts. 

In accordance with a fifth aspect of the invention, there is provided a wavefront 
sensor comprising: 
imaging means 
10 aberration means 
detecting means 
feedback means, 

arranged such that the imaging and detecting means capture images of the wavefront 
field in a manner such that the wavefront field imaged is subject to two or more 
15 different aberrations imposed by the aberration means. 

Preferably, the aberration means when combined with the imaging means can deliver 
to the detecting means simultaneously or sequentially images of the input wavefront 
field subject to different but known aberration point spread functions (psf) such that 
20 the filter functions (Fourier transforms of the aberration psfs) have complex odd or 
complex even symmetry and at least two of the images are subject to filter functions 
that are complex conjugate pairs. 

The aberration means does not provide pure defocus aberrations. 

25 

The filter functions may be described by any pure Zernike polynomial. 

Preferably, an electronic feedback means is arranged to deliver to an adaptive optics 
wavefront modulator a signal derived from the difference between the images affected 
30 by the complex conjugate filter functions and used to drive the wavefront modulator 
in such a manner that the feedback signal strength is zero when the images recorded 
with complex conjugate filter functions are identical. 



WO 2004/113856 PCT/GB2004/002657 

9 

Preferably, the filter functions can be a weighted sum of Zernike polynomials 
arranged optimally to equalise the signal generated from each mode of deformation to 
be corrected in the input wavefront field according to the expected statistical 
distribution of such modes in the input wavefront field. 

Preferably, the aberration means is a diffractive optical element arranged such that the 
complex conjugate filter functions are associated with diffraction orders of the same 
order but different sign (±1, ±2, etc). Alternatively, the aberrating means can be a 
variable-shape optical mirror or a variable refractive index device such as a liquid 
crystal phase modulator used sequentially to provide complex conjugate filter 
functions. Optionally, the aberration means can also consist of a deformed reflective 
surface where the illumination of that surface from each side produces the complex 
conjugate filter functions. 

The following analysis show that the filter functions created by the aberration means 
are capable of encoding the sense and location of wavefront errors on heavily 
scintillated and discontinuous input wavefronts. The necessary and sufficient 
conditions required of filter functions suitable for use in a generalised phase-diversity 
null wavefront sensor for adaptive optics applications are quantified. 

Let Y(r) = | x P (r)| e l<p(r) represent the complex amplitude distribution in the entrance 
pupil of an optical system, r being the co-ordinate in the pupil plane. If ^(r) 
represents a plane wavefront, the phase satisfies <p (r) = constant . Any spatial 

variation in the wavefront phase represents a distortion, or aberration, that requires 
correction in an AO system. 

To operate successfully as a null wavefront sensor in an AO system we require a 
device that produces an error signal where, <p (r ) is not constant. For the purposes of 
this analysis, the optical output from monochromatic systems are considered. Thus a 
discontinuity of exactly an integral number of wavelengths in size does not affect the 
optical output, e.g. image intensity profile, and may be regarded as an undistorted 
wavefront. Such a discontinuity may be an important wavefront error in metrology 
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applications. It is desirable that the error signal provides information that localises the 
wavefront error in r-space and indicates the sense in which correction should be 
effected (in the absence of such indication a multi-dither technique™ 11 is required to 
effect correction). 

5 

Let ^(£)=i/(£) + ,4(£) be the Fourier transform of ¥(r), where H(4) and A(%) 
represent respectively the Fourier transforms of the real and imaginary parts of ^(r). 
Clearly, H(g) is Hermitian and A(£) is anti-Hermitian. These symmetry properties 
will be required later. Thus, 

10 KhJ ; h) (1) 

Let F ± (£) « R (£) ± H be the filter function that is the Fourier transform of an 
aberration psf f ± (r) , which represents a complex function with which *F (r) is 
convolved when forming an image of the system entrance pupil. The functions R (£) 
15 and are real-valued functions and the± indicates the use of two filter functions, 
in which the Fourier phase of the filter is reversed. In this work we are particularly 
interested in the necessary and sufficient conditions that constrain R(4) and in 

such a way that F ± (£) are suitable filter functions to provide a null wavefront sensor 
for use in adaptive optics. 

20 

The physical relationship between these different functions is indicated schematically 
in figure 2 where wavefront 1 1 is shown along with a diffractive optical element 
which provides the aberration means 13. Images 15 and 17 with aberration function j 
+(r) and functions j.(r) respectively are shown along with an image 19 of the 
25 wavefront. 



The detected intensity function may thus be written 



(2) 
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Substituting for F ± , expanding and simplifying, the difference between the images 
formed using the two filter functions may be expressed 
d(r) = j + (r)-j_(r) 

5 This is a real- valued function, since the quantity in [ ] is a difference of two complex 
conjugates and is thus imaginary-valued, so 

21 (4) 

The rhs of (4) can then be expanded and the terms grouped into 4 separate expressions 
10 which are equal to the rhs of (4) when summed: 

J^jy(^)/(^)e^ J^'iT (^')^(^ - J^^)^)^ J^'^(^ ')/(^')^' (5.1) 

f rf^^c^)^^)^ ^ J C^O^C^O^ ^ - J rf^^4C^>^(^>^"^ J rf^'^* C^O^C^O^" ^ (5.2) 

15 J^(#)J(£K*^ (5.3) 



This expression for the difference between the two detected intensity functions is 
20 generally valid - no restricting assumptions have so far been made. We may now 
investigate under what conditions of symmetry these expressions individually, or 
summed, are identically zero. 

Conditions for operation as a null wavefront sensor 

25 

3 .2. 1 Filter function must be complex 
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Unless both R and / are non-zero, all of the terms in equations (5) are identically 
zero V \jf . Thus d (r) is identically zero for all input wavefronts and no error signal is 
generated from a non-flat wavefront. 

5 This may be readily understood since if I(Q is zero F ± is identical whatever the 
arithmetic sign, so the two images are identical Vys . If R(Q is zero the change of 
arithmetic sign is lost by the modulus square operation that is an inevitable part of 
quadrature detection processes - and thus the two images are again identical V y/ . 

10 3.2.2 Complex filter function 
3.2.2.1 Even symmetry 

Suppose that both and R{£) are even functions of £ . 

15 

Consider the first expression (5.1). Since H{g) is Hermitian and is symmetric 
and real-valued the product H(%)I(%) is Hermitian. Thus the Fourier integral 
fd£HffiI(g)e-* 4 is real-valued. The same is true of $d£'H\Z')R(t;')e ir * . The 

second product of two integrals is term by term the complex conjugate of the first 
20 product. Thus expression (5.1), the difference between two complex conjugates, is 
always zero when both and are symmetric. 

Similarly (5.4) is always zero because each of the integrals reduces to a purely 
imaginary function. The product of these imaginary functions is real and the 
25 difference between the two complex conjugate terms is again always zero. 

Thus the difference between the two images is the sum of (5.2) and (5.3) and can be 
written 
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(6) 

In equation (6), for each pair of integrals one integral reduces to a real-valued function 
5 and the other integral to an imaginary-valued function. The first and second terms, 
also the third and fourth terms, are complex conjugate pairs and thus the rhs of (6) is 
imaginary valued or zero. 

If either H or A is zero or if A = He l(p with <p constant, equation (6) and thus d (r) , 

10 is zero. However, these conditions are exactly those under which the input wavefront 
is flat and a null sensor is required to produce a null output. 

Thus, if the filter function is complex with even symmetry the difference between the 
two images formed using these filters forms a potentially-useful null wavefront 
15 sensor. 

3.2.2.2 Odd symmetry 

Suppose that both I (<f ) and R (£) are odd functions of £ . 

20 

In (5.1) and in (5.4) the odd symmetry of the real functions R and I means that the 
products within the integrals in (5.1) have anti-Hermitian symmetry and those within 
(5.4) have Hermitian symmetry. Thus, although the arguments given are reversed for 
each term from those given in 3.2.1, both (5.1) and (5.4) are identically zero Vy/ . 

25 Thus d (r) again reduces to equation (6). The arguments from 3.2.1 again hold, 

although the role of the terms is reversed, one term in each integral product is purely 
real and the other term purely imaginary. 

Thus a filter function with complex odd symmetry is potentially suitable for use as a 
30 filter function for wavefront sensing. 



WO 2004/113856 

3.2.2.3 Mixed symmetry 



14 



PCT/GB2004/002657 



Suppose that one of / (£) and R (£) is an even function of £ and the other is an odd 
functions of £ . 

5 

In (5.1) the mixed symmetry will result in one of the integrals in each product being 
purely imaginary and the other integral purely real (dependent on whether I or R is 
odd). In either case the product of the integrals is purely imaginary and thus (5.1) is 
purely imaginary or zero. An equivalent argument shows that (5.4) is purely 
10 imaginary or zero. 

Note, however, that the expressions (5.1) and (5.4) are reliant on the interactions 
between H and the filter function or on A and the filter function and do not involve 
cross terms between H and A . 

15 

In equation (6) the mixture of odd and even symmetry will mean that both integrals in 
each product are either purely real or purely imaginary. In each case the product of 
these terms will be purely real and the difference of the complex conjugates will thus 
be zero Vy. 

20 

Thus, filters with mixed symmetry produce no error signal dependent on the deviation 
of the wavefront from a plane wave and thus such filters are unsuitable for use as a 
wavefront sensor. 

25 3.2.3 Sensing the error direction 

Information about the wavefront aberrations is encoded by the sum of the cross terms 
in equation (6). 

30 If the sense of the wavefront error is reversed, the phase of the wavefront will change 
arithmetic sign. This means that for a wavefront error of a given amplitude the error 
signal (equation (6) for filter functions with complex even or complex odd symmetry) 
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changes sign if the sense of the error is reversed. The wavefront sensor thus delivers 
an error signal that preserves information about the sense of the wavefront error. 

Note that, since the relationship between the error signal and the wavefront error is 
5 non-linear (involving the balance between H and A in the description of the 

wavefront), this does not guarantee that the error signal can be inverted to find the 
wavefront shape. 

3.2.4 Localising the wavefront error 

10 

Returning to equation (6) we note that each of the integrals, when expressed as a 
function of r appears in the form of a convolution of H or A with a real-valued 
function deriving from the filter function. As the products are of complex functions 
and real-valued function the phase of the product is identically equal to that of the 
15 Hermitian or Anti-Hermitian function in each integral. 

Without loss of generality, the location of the wavefront error can be identified with 
the position at which a(r) , the Fourier transform of A(£) , is non-zero. Because no 

phase ramp is added by the product with the real-valued filter function, the location of 
20 the non-zero component is, in some sense, localised at the point where a(r) is non- 
zero. 



The key term that encodes information about the wavefront aberrations is the sum of 
the cross terms in equation (6). 

25 

If the sense of the wavefront error is reversed the arithmetic sign of A (£) reverses 
because the phase of the wavefront will change arithmetic sign. This means that, for a 
wavefront error of given amplitude, the error signal (equation (6) for filter functions 
with complex even or complex odd symmetry) changes sign if the sense of the error is 
30 reversed. The wavefront sensor thus delivers an error signal that preserves 
information about the sense of the wavefront sensor. 
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Equation (6) also provides localisation of the wavefront error. Each of the integrals, 
when expressed as a function of r, appears in the form of a convolution of the Fourier 
transform of H or A with a Hermitian function related to the filter function. Without 
loss of generality, the location of the wavefront error can be identified with the 
5 position at which a(r), the Fourier transform of A{£) , is non-zero (a(r) is directly 

related to the phase distortion). The function with which the convolution takes place 
may or may not have a maximum at the origin. When it does not the location of the 
non-zero component is, in some sense, localised at the point where a(r) is non-zero, in 
the sense that it may be distributed around the position at which a(r) is non-zero. 
10 When the intensities are recorded in pupil space this provides the localization required 
for a null wavefront sensor. 



The present invention will now be described by way of example only with reference 
to the accompanying drawings of which: 

15 

Figures la and lb show a series of wavefronts propagating and illustrates variations in 
local beam intensity; 

Figure 2 illustrates a wavefront and a diffractive optical element which provides the 
aberration means; 
20 Figure 3 illustrates an embodiment of the present invention; 

Figures 4a to 4d show simulation result with the non-plane wavefront input and the 
even symmetry phase diversity filter; 

Figure 5 shows a second embodiment of the present invention; 
Figure 6 shows a third embodiment of the present invention; 
25 Figure 7 shows a fourth embodiment of the present invention; 
Figure 8 shows a fifth embodiment of the present invention; 

Figure 9 shows results with the Zernike polynomial Zio° phase diversity filter used to 
create aberration means; 

Figures 10(a) (i) to (iii) show the wavefront intensity of a wavefront in the pupil 
30 phase, the output intensity distribution of a difference image and the associated prior 
art quadratic diffractive optical element; 

Figures 1 1(b) (i) and (iii) show an embodiment of the present invention and illustrate 
the wavefront intensity of a wavefront in the pupil phase, the output intensity 
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distribution of a difference image and the associated pure spherical diffractive optical 
element; and 

Figures 1(a) to 1 1(b) show a discontinuous pupil phase wavefront and outputs with a 
prior art defocus grating and a grating according to the present invention. 

5 

The requisite filter functions could be implemented in several ways, but one of the 
most straightforward approaches (for use with quasi-monochromatic radiation) is to 
use a diffractive optical element (DOE) to introduce the aberration. 

10 Figure 2 shows a DOE used to implement the requisite filter functions in the + 
diffraction orders and the zero order contains an image of the field to be 
reconstructed. In situations where the photometric efficiency is important a phase 
grating may be used and the depth of the grating adjusted to reduce to zero, or close to 
zero, the flux in the zero diffraction order. In Figure 2, the DOE 13 is shown that is 

15 used to generate the aberrations, the image with aberration psf j+(r), 15 and the image 
with aberration psf j.(r), 17 are also shown along with a representation of the 
wavefront to be reconstructed 19. From 19 we can tell whether the wavefront to be 
reconstructed is scintillated, multiply connected, uniform intensity etc. 

20 Figure 3 shows the combination of DOE 27 (diffraction grating) with SLM 29 (spatial 
line modulator) operated using liquid crystals in this example. The wavefront 
modulators are used to produce an adaptive optical system for driving a null sensor. 
In Figure 3 the CMOS (complementary metal oxide semiconductor) camera 23 used 
to detect the j+1 and j-1 signals is shown along with the lens 27and logic circuit 31 

25 which provides an electrical control system that is used to control the SLM 29 thereby 
altering the shape of the wavefront in response to errors detected in the j+1 and j-1 
signals. 

In this example, the difference between the two images on a pixel by pixel basis is 
30 used directly to drive the wavefront modulators in order to provide the requisite 
wavefront correction. 



Spatial Light Modulators (SLM's) can be used to provide the wavefront modulation. 
SLM's provide a low-cost, low-voltage and low-volume modulator. In addition, 
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diversity can be applied to the input wavefront in the grating and lens combination. 
The grating is provided with a carefully selected filter function that satisfies the 
necessary and sufficient conditions for a null wavefront sensor in accordance with the 
present invention. The grating allows the source image (in the 0th order) and the 
5 phase diverse data (in the ±1 orders) to be displayed on the focal plane, in which a 
CMOS camera is located. 



The images in the diffraction orders are subtracted to deduce the error signal and this 
information is fed back to the modulators to adjust the correction of the incoming 
10 wavefront. 

A number of examples of the effect of a null wavefront sensor in accordance with the 
present invention are now provided. The above analysis has shown that a wavefront 
sensor can be created using a range of aberration means typically formed from 

1 5 diffractive optical elements. Furthermore, phase diversity filters can be used to give 
superior performance to the defocus method by using a priori information about the 
wavefront to be sensed. As seen in figure 2, image 19 shows the type of wavefront 
where the wavefront is scintillated, discontinuous or has disconnected wavefront 
segments. Aberration means is accordance with the present invention are used in the 

20 wavefront sensor to provide an accurate tool for determining wavefront error and 
therefore correcting the error in the wavefront. 

The embodiments of the present invention shown in Figures 3, and 5 to 8 of the 
present invention can be described as Compact adaptive optical systems (CAOS). In 

25 general, the apparatus is based on transparent wavefront modulators (e.g. liquid 

crystal SLMs) and phase diversity wavefront sensing. As described above, Compact 
Adaptive Optical (AO) systems, based on the use of DOEs, phase diverse wavefront 
sensing and transparent wavefront modulators, offer the potential for construction of 
AO systems with benefits such as: 

30 1 • Minimisation of non common-path errors by combining the wavefront 

sensor data and the corrected image in a single focal plane with essentially 
no separation of the 'science' and wavefront sensing optical trains; 
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2. The ability to place wavefront modulators conjugate to multiple planes in 
the object space whilst preserving a compact and robust optical train that 
would be particularly well-suited to multi-conjugate; and 

3. Avoidance of the requirement for conventional optics, such as beam 
5 splitters and conventional bulk optical components for re-imaging 

conjugate planes, that increases the size and weight of the equipment. 

Experimental simulations have been undertaken with plane and non-plane input 
wavefronts and a range of phase diversity filter functions having properties in 
10 accordance with the present invention. 

Figures 4a to 4d show simulation result with the non-plane wavefront input and the 
even symmetry phase diversity filter (Zi 0 °, Zemike polynomial). From the simulation 
result we can see that the intensity difference between +1 order and -1 diffraction 
15 images is non-zero for the non-plane wavefront input. This is a sufficient condition 
for a null wavefront sensor derived from the theory. Figures 4 correspond to data 
recorded in an experimental embodiment of the form shown in figure 8. 

Figures 5 to 8 show additional embodiments of the present invention. 
20 For reasons of clarity, the detector/camera has been omitted from these figures. 

Figure 5 shows an embodiment 35 of the present invention in which an input/ test 
wavefront 37 is shown at a position remote from a lens/grating 39. A modulator 41 is 
also shown connected to the output 51 of the logic circuit 49 which is used to 
25 calculate the difference in the intensity signals generated at l x and I 2 denoted by 

reference numerals 43 and 45 respectively. Io describes the position of the image of 
the wavefront 37. 

In this example of the present invention, the input/test wavefront is positioned 
30 remotely from the lens/grating 39. It has been found that such an arrangement 

provides a sufficient degree of spatial resolution of the Ii and I 2 signals to allow the 
actual position in the wave front of the non-planar part to be identified. This feature is 
particularly useful where the present invention is to be used as an adaptive optical 
device for the correction of input wavefronts . In this example, the measured , 
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spatially separated values of Ii and I 2 will contain dark and bright areas where the 
wavefront has deviated from being a planar wavefront. Once the difference between 
the calculated values of Ii and I 2 has been calculated, the logic circuit 40 will create a 
signal which allows the modulator to adjust the shape of the wavefront to make it 
5 planar. 

Figure 6 shows a further embodiment 53 of the present invention in which the 
input/test wavefront 55 is arranged in the plane of the lens/grating 39. In this 
example, the output at Ij and I 2 provide less well spatially resolved signal. Therefore, 
it is more difficult to obtain information on the exact position across the wavefront 
10 where the wavefront is non-planar. In this example the logic circuit is used to provide 
a good indication of the size of the error, but less spatial information. A value for the 
overall wavefront error may be obtained relatively quickly. Figures 4 and 9 show 
simulations of the data obtained for embodiments of the type shown in figure 6. 

15 Figures 7 and 8 provide additional embodiments of the present invention with figure 7 
showing a wavefront remote from the lens/grating and figure 8 showing a wavefront 
in the plane of the lens/grating. In addition, a modulator 79 is provided in the 
apparatus of figure 8 to allow resolution of complex wavefront shapes. 

20 The examples of figures 7 and 8 are primarily intended for use with analytical 

instruments where the wavefront shape 59, 75 is a reflection from a surface such as a 
mirror. The reflections are used to measure the shape and quality of the surface from 
which the wavefront has been reflected. As with the arrangement shown in figure 5, 
figure 7 provides an output from Ii and I 2 that is relatively highly spatially resolved 

25 and which provides an output 51 in the form of data on the position of the wavefront 
error, which in turn provides information on the defect in the surface. As with the 
arrangement in figure 6, the arrangement in figure 8 provides a relatively quick 
measure of the magnitude of an error. 

30 Typically, these instruments can be used for measuring defects in astronomical 
mirrors. 
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Figure 9 shows results with the Zernike polynomial Zio° phase diversity filter 
function used to create aberration means in an embodiment of the present invention. 
This filter function is an even phase diversity. 

5 Figures 9a and 9b show the original signal phase and the retrieved signal phase. 
Figure 9c shows the difference between these phases, figure 9c shows the In error 
value plotted against the number of iterations and figures 9e and 9f show the 1+1 and 
1-1 images. 

10 Figures 1 0(a) (i) to (iii) show the wavefront intensity of a wavefront in the pupil 

phase, the output intensity distribution of a difference image and the associated prior 
art quadratic diffractive optical element. Figures 10(b) (i) to (iii) show the wavefront 
intensity of a wavefront in the pupil phase, the output intensity distribution of a 
difference image and the associated pure spherical diffractive optical element, an 

15 embodiment of the present invention. The prior art defocus grating is constructed 
using a curve described by the expression 2^ 2 - 1 . The spherical grating is 
constructed using a curve described by the expression 6£ 4 - 6£ 2 + 1 . 

The mathematical shape of the input wavefront has sufficient continuity and 
20 uniformity of intensity to operate with the prior art defocus grating and a grating 

suitable for use in the present invention.. As shown in the diagram of the pupil phase 
in Figures 10(a)(i) and 10(b)(i) are substantially identical. However, it will be noted 
that the amplitude of the input wavefront error for figure 10(b)(i) is approximately 
10% of the amplitude error shown in figure 10(a)(i). The difference images 10(a)(ii) 
25 10(b)(ii) have substantially, the same intensity. Therefore, the wavefront sensor 

created using filter functions which are complex valued and which have non-mixed 
symmetry has superior sensitivity compared to the defocus only quadratic filter. 

Other non-quadratic polynomials maybe used to define curves for constructing 
30 gratings in accordance with the present invention. 

The error graph is created using a known iteration algorithm and indicates that the 
error reduces to a substantially constant value after around 50 iterations in this 
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example. It has been found that the error values for the Zio° polynomial are in this 
example, lower that those achieved using a defocus filter. 

Figures 11a and 1 lb show wavefront intensity in the pupil phase and output intensity 
5 after the wavefront has passed through a defocus gravity as shown in figure 10(a)(iii). 
Figures 11c and lid show wavefront intensity in the pupil phase and output intensity 
after the wavefront has passed through a spherical aberration diffraction gravity as 
shown in figure 10(b)(iii). In both figures 1 1(a) and 1 l(c) 5 the input wavefronts are 
identical and have a central obscuration 50 in the centre of the wavefront. The central 
10 obscuration shown is of the type found in telescopes. 

It can be seen from figure 1 1(b) that there are a number of features around the 
boundary of the central obscuration 52 when the wavefront passes through a defocus 
gravity. These features arise because of the sensitivity of the defocus gravity to the 
15 missing boundary and do not represent errors in the ^avefront, but are created by the 
optical properties of the defocus gravity. In a wavefront sensor, the presence of such 
features would prevent accurate modification of the wavefront shape and prevent the 
creation of a plane wavefront. 

20 It can be seen from figure 1 1(d) that the boundary around the central obscuration 54 is 
smooth. Therefore, the use of a grating in accordance with the present invention is 
which the grating is constructed such that the shape of the aberrating means is defined 
by a filter function that is complex valued and has non-mixed symmetry. 

25 In further examples of the present invention, the input wave can be moved and the a 
grating whose shape is defined by a polynomial in accordance with the present 
invention, selected to provide optimum results. This can be achieved through 
statistical analysis of the wavefront. 

30 Improvements and modifications may be incorporated herein without deviating from 
the scope of the invention. 
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